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Abstract. The three-dimensional evolution of a pure electron plasma is studied by means of a 
particle-in-cell code which solves the drift-Poisson system where kinetic effects in the motion 
parallel to the magnetic field are taken into account. Different results relevant to the nonlinear 
dynamics of trapped plasmas and low-energy electron beams are presented. 

INTRODUCTION 

The evolution of an electron plasma in a Malmberg-Penning trap [1] is studied by means 
of the particle-in-cell (PIC) code MEP [2]. The electron dynamics is analyzed in the 
frame of a guiding center electrostatic approximation, where the velocity perpendicular 
to an externally applied uniform axial magnetic field is given by the electric drift, and 
kinetic effects in the motion parallel to the magnetic field are taken into account. The 
evolution of the system is followed within a conducting cylindrical surface of radius 
R and length L on which the (in general time-dependent) boundary conditions for the 
electrostatic potential are imposed. 

The code is applied here to investigate two different situations (see Fig. 1). In the first 
case, the evolution of a traveling electron beam is considered. It is observed in particular 
how in a space-charge dominated regime a fraction of the electrons close to the axis 
is reflected back to the cathode while a high density annulus is formed inside the drift 
tube, which gives rise to the development of vortex-like structures. In the second case, 
the injection phase of an electron plasma in a Malmberg-Penning trap is studied. Here, 
the electrons enter in an initially empty trapping region, are reflected by a fixed potential 
barrier on the opposite side of the trap and come back to the cathode, thus interacting 
with new emitted electrons. It is shown how, in dependence of injected current and 
geometrical parameters, a virtual cathode may form close to the injection surface, and 
a hollow electron column is formed in the trapping region. Longitudinal kinetic effects 
are also investigated by varying the velocity distributions of the injected electrons. 



MODEL AND PIC CODE 

The system is described in the frame of the zeroth order drift approximation, where the 
velocity perpendicular to the uniform magnetic field B = Be z (e z being the unit vector in 
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FIGURE 1. Scheme of traveling beam configuration (left) and trap filling configuration (right). —V is 
the trap potential. 



the z direction) is given by the electric drift, \e = (c/B) e z x V<p, with <p the electrostatic 
potential and c the speed of light. Assuming the guiding centers distribution, /, of the 
form f(x,\,t) = F(v|i,x,f) 8(v± — \e), where x,\,t are position, velocity and time, 
respectively, 8 denotes the Dirac's distribution, and vn and Vj_ denote the components of 
the velocity parallel and perpendicular to the magnetic field, the Vlasov-Poisson system 
reads 

dF 1 1 dF 

+ ( e z xV<p +V ||e z )-VF + -— e z -V<p — = ; V 2 <p = n, (1) 

at 2 11 M e ff ovii 

where n(x,t) = J F dv\\ is the electron density. Adimensional quantities are used: length, 
time, density and potential are normalized over R, 0) c /2ft>~, hq and AnenoR 2 , respec- 
tively, where G) c = eB/mc is the non-relativistic electron cyclotron frequency, co p = 
(4-7te 2 no/m) 1 / 2 is the electron plasma frequency, computed for a specified electron den- 
sity no, and —e, m are electron charge and mass, respectively. The behavior of the system 
is therefore characterized by the single parameter M e /f, which plays the role of an ef- 
fective mass, 

M eff = 4^ 4.115- (2) 
n co 2 5 2 [kGauss] 

It results M e ff = Iuq/ub, ng being the so-called Brillouin density, n# = (B 2 /Sn)/mc 2 , so 
that < M e ff < 2. Introducing a cylindrical system of coordinates (r, 6,z), and using the 
variable s = r 2 , the equations of motion corresponding to the trajectories of the kinetic 
equation in (1) are: 

ds _ d(p dQ _d(p dz _ dv\\ 1 d<p 
dt dd ' dt ds ' dt " ' dt M e /f dz ' 

In the code, Eqs. (1) are discretized on an equispaced grid (with the only exception 
of the central cell) in the coordinates s, 6 and z. The number of cells is denoted as 
N s , Nq and N z , respectively. The grid for s is defined as s_ = 0, = 1/(N s Nq + 1), 
= lj+N e /(N s Ne + 1), so = 0, sj = (sj + s j+l )/2, j = l...N s (s is the position 
of the center, while s denotes the lower boundary of a "radial" cell). The grid for G 
is Oj = 2n(l — 1)/Nq, 1 = 1 ...Nq, with the periodicity relation 0n 8 +i = 9i, while the 
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FIGURE 2. Left: Time evolution of a beam in the plane (z, r). The parameters of the run are: M e f f = 
0.01, N s = 64, N e = 64, N z = 128, N p = 10 s . The beam is injected at z = -10, with an initial parallel 
velocity viiq = 9.0. The current is / = 0.2215, / being normalized over the ratio between enR 2 Lno and the 
normalization time. Emitter and collector are both at ground. From left to right, top to bottom, the data 
refer to t = 0.7, 1 .5, 3.5 and 10.0, respectively. Right: Transverse projection [(x,y) plane] of the particles 
in the interval — 10 < z < —5, — 5 < z < 0, < z < 5 and 5 < z < 10, respectively. The data refer to 
f = 10.0. 



grid for z is Zk = (k — 1 /2)L/N Z — L/2, k = 1 ...N z . Each cell has the same volume 
AV = AsAOAz/2, with As = l/(N s + 1 /N e ), AO = 2n/N e , and Az = L/N z , respectively. 
The system governed by Eqs. (1) is simulated numerically as an ensemble of macro- 
particles with fixed sizes As, AO and Az, using a PIC method [2]. 



RESULTS 

The code is used to simulate the plasma dynamics in a Malmberg-Penning trap for 
different possible experimental settings. The physical situation is determined by several 
geometrical and emission parameters: the magnetic field strength B; the sizes R, L and 
the geometry of the emitting surface; the potentials which are imposed on cathode, anode 
and drift tube; the initial electron velocity distribution and the initial current distribution 
emitted by the source. Malmberg-Penning traps usually use a spiral-wound tungsten 
filament for the injection [3]. The MEP code is able to simulate this initial spatial 
distribution of the electrons; the effect of an accelerating grid is considered by suitably 
"cutting" the spiral along rows and columns of a given width. In addition, the code is 
able to take into account various initial velocity and current density distributions. At 
first, a low-energy electron beam traveling in an equipotential (grounded) drift tube is 
considered (see Fig. 1 left). The electrons are continuously emitted from the cathode and 
collected to the anode. The characteristic time scale of the system is the time of flight of 
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FIGURE 3. Left: Time evolution of a plasma filling the trap, in the (z,r) plane. The parameters of the 
run are: M eff = 4.64 • 1(T 4 , N s = 64, N e = 64, N z = 128, N p = 10 5 . The beam is injected at z = -5.0, 
with V|| = 28.90. The current is / = 0.1423, and the potential difference between collector and emitter is 
—0.97. From left to right, top to bottom, the data refer to t = 0.7, 1.5, 4.0 and 8.0, respectively. Right: 
Transverse projection of the particles in the interval — 5 < z < —2.5, —2.5 <z<0, 0<z< 2.5 and 
2.5 < z < 5, respectively. The data refer to t = 8.0. 



an electron, L/v\\. In a space-charge-limited regime, it is found experimentally that the 
central part of the beam is reflected, a hollow electron column forms and fast coherent 
structures arise, possibly due to the development of diocotron instability [4]. The PIC 
simulations confirm the experimental findings. The time evolution of a mono-energetic 
beam in the (r, z) plane is shown in Fig. 2 left. It is evident, in particular, the formation of 
a virtual cathode close to the injection surface. The central part of the beam is reflected 
back to the cathode by the space-charge of the beam itself, and only the outer part of 
the beam reaches the opposite end of the trap. This annular beam shows a quasi-2D 
evolution of vortex structures: Fig. 2 right represents the electron plasma distribution in 
(almost) stationary conditions on different transverse slices. 

As a second example, the case of an electron plasma filling the trap is studied (see 
Fig. 1 right). This situation simulates the phase of injection in a Malmberg-Penning trap. 
In this phase, the end plug electrode is maintained at a sufficiently negative potential in 
order to reflect the electrons entering the trap, while the plug electrode at the entrance 
is grounded to let the electrons flow into the trap from the cathode (starting from this 
configuration, the trapping phase is obtained by simply lowering the plug potential at 
the entrance of the beam to the same value as the end potential, so that the electrons 
are electrostatically confined along the axis of the device). A mono-energetic beam is 
considered. It is found that when the input current of the beam is increased, the space 
charge cloud which forms inside the trap represents a barrier not only for the new 
incoming particles, but also for the electrons which are already inside and are reflected 
by the external negative potential barrier. As a result, the electron plasma filling the trap 
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FIGURE 4. Time evolution of a plasma filling the trap, in the (z,p\\) plane. The parameters are 
the same as in Fig. 3 (red dots). From left to right, top to bottom, the data refer to t = 0.3, 0.7, 
1.5 and 3.0, respectively. The blue dots are relevant to an initial Maxwellian distribution in parallel 
velocity, F(vii,x,0) = g(x) -exp[— (v|| — V|| ) 2 /2(7 2 ), with the same central value V||q, and a velocity spread 
(7 = 5.78. 



assumes an annular shape, as shown in Fig. 3. The time evolution of the electrons in 
the phase plane (z,p\\) is shown in Fig. 4 (red dots), where pn = M e ffVn. Longitudinal 
kinetic effects have been investigated in this case. In particular, it has been found that 
varying the parallel velocity distribution of the injected electrons to a Maxwellian has 
a dramatic effect (see again Fig. 4, blue dots): space-charge effects are much weaker in 
this case, and the electron plasma maintains its initial transverse (spiral) shape. 
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